FERMILAB-PUB-09-102-A-T 
UK/TP 09-02 

Dark Matter Constraints from a Cosmic Index of Refraction 

Susan Gardner^'^* and David C. Latimer^ 

^Center for Particle Astrophysics and Theoretical Physics Department, 
Fermi National Accelerator Laboratory, Batavia, IL 60510 and 

'^Department of Physics and Astronomy, 
University of Kentucky, Lexington, KY 40506-0055 

Abstract 

The dark-matter candidates of particle physics invariably possess electromagnetic interactions, if 
only via quantum fluctuations. Taken en masse, dark matter can thus engender an index of refrac- 
tion which deviates from its vacuum value. Its presence is signaled through frequency-dependent 
effects in the propagation and attenuation of light. We discuss theoretical constraints on the ex- 
pansion of the index of refraction with frequency, the physical interpretation of the terms, and 
the particular observations needed to isolate its coefficients. This, with the advent of new oppor- 
tunities to view gamma-ray bursts at cosmological distance scales, gives us a new probe of dark 
matter and a new possibility for its direct detection. As a first application we use the time delay 
determined from radio afterglow observations of distant gamma-ray bursts to realize a direct limit 
on the electric-charge-to-mass ratio of dark matter of |e|/M < 1 x 10^^ eV^^ at 95% CL. 
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Some twenty-three percent of the Universe's energy budget is in dark matter [l|-|3|], yet, 
despite its abundance, httle is known of its properties. A number of methods have been 
proposed for the detection of dark matter. Such studies typically rely, for direct searches, 
on dark-matter-nucleus scattering, and, for indirect searches, on two-body annihilation of 
dark-matter to Standard Model (SM) particles; constraints follow from the nonobservation 
of the aftermath of particular two-body interactions. In contrast, we probe dark matter in 
bulk to infer constraints on its particulate nature. 

We search for dark matter by studying the modification of the properties of light upon 
passage through it. One can study either polarization j3] or propagation effects; we focus 
here on the latter. The resulting constraints are most stringent if dark matter consists of 
sufficiently low mass particles, be they, e.g., warm thermal relics or axion-like particles, that 
its number density greatly exceeds that of ordinary matter. We thus consider dark matter. 
Matter effects are signaled by dispersive effects in the speed or attenuation of light. We study 
this by introducing an index of refraction n{u,z), whose deviation from unity is controlled 
by the light-dark- matter scattering amplitude in the forward direction, i.e., the forward 
Compton amplitude, as well as by the angular frequency cu of the light and the redshift z 
at which the matter is located. A dark-matter particle need not have an electric charge 
to scatter a photon; it need only couple to virtual electromagnetically charged particles to 
which the photons can couple. The scattering amplitude is related by crossing symmetry 
to the amplitude for dark-matter annihilation into two photons, so that any dark-matter 
model which gives rise to an indirect detection signal in the two-photon final-state [5|] can 
also drive the index of refraction of light from unity. Its real part is associated with the 
speed of propagation, and we search for its deviation from unity by searching for frequency- 
dependent time lags in the arrival of pulses from distant gamma-ray bursts (GRBs). 

The limits on the non-observation of frequency-dependent effects in the speed of light 
are severe. The best limits from terrestrial experiments control the variation in the speed 
of light c with frequency to |5c|/c < 1 x 10~^ j6|, but the astrophysical limits are much 
stronger. The arrival time difference of pulses from the Crab nebula bound |5c|/c < 5 x 
10"^'' 0, Isj, and still better limits come from the study of GRBs jl, [oj. GRBs are bright, 
violent bursts of high-energy photons lasting on the order of thousandths to hundreds of 
seconds, and their brightness makes them visible at cosmological distances. Time delays 
which are linear in the photon energy can occur in vacuo in theories of quantum gravity; the 
special features of GRBs make them particularly well-suited to searches for such signatures 
of Lorentz violation [9]. The detection of photons of up to ~ 31 GeV in energy from GRB 
090510 severely constrains this scenario, placing a lower limit on the energy scale at which 
such linear energy dependence occurs to 1.2 Epianck ~ 1-5 x 10^^ GeV [lo|. Although vacuum 
Lorentz violation and light-dark-matter interactions can each induce dispersive effects, their 
differing red shift and frequency dependence render them distinct. 

A model-independent analysis of the deviation of the refractive index from unity is possi- 
ble if we assume that the photon energy is small compared to the energy threshold required 
to materialize the electromagnetically charged particles to which the dark matter can cou- 
ple. In models of electroweak-symmetry breaking which address the hierarchy problem, a 
dark- matter candidate can emerge as a by-product. In such models the inelastic threshold 
Uth is commensurate with the weak scale, or crudely with energies in excess of (9(200 GeV), 
as the new particles are produced in pairs. If the photon energy u satisfies the condition 
uj <^ Uth, we can apply the techniques of low-energy physics to the analysis of the forward 
Compton amplitude. Under an assumption of Lorentz invariance and other symmetries, we 
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expand the forward Compton amplitude in powers of oj and give a physical interpretation 



to the coefficients of the first few terms as w — >■ llMlSl]. In particular, the term in 0{u:^) 
is controlled by the dark matter particles' charge and mass, weighted by preponderance, 
irrespective of all other considerations save our assumption of Lorentz invariance • 

The relationship between the index of refraction njuj ) an d the forward scattering ampli- 
tude /cj(0) for light of angular frequency to is well-known [l5|,ll6|, where we relate /^^(O) to the 
matrix element Ai of quantum field theory to connect to particle physics models of dark mat- 
ter. Using standard conventions [13], we determine n{u:) = 1 + {p/4M'^u'^)Air{k,p k,p) , 
in the matter rest frame [l^, so that p = (M, 0) and k = {u, uh) with p the mass density of 
the scatterers and M the particle mass. In our analysis we assume M ^ T{z), where T{z) 
is the temperature of the dark matter at the red shift of the observed gamma-ray burst. 
Since T{z) should be a factor of some ((1 + z)/{l + 2;prod))^~^ smaller than T(2;pi.od) at the 
moment of its production or decoupling, our limits are not restricted to cold dark matter ex- 
clusively. Moreover, even in the latter case, the candidate mass can be as light, e .g., as light 



as M ~ 6 X 10~® eV in the axion model with Bose-Einstein condensation of Ref. Under 
the assumptions of causality or, more strictly, of Lorentz invariance, as well as of charge- 
conjugation, parity, and time-reversal symmetry in the photon-dark-matter interaction, we 
have [Ul, [l3| A4r{k,p k,p) = fi{uj)e'* ■ e + if 2(00)8 • e'* x e, where S is the spin operator 
associated with the dark-matter particle and e (e') is the polarization vector associated with 
the photon in its initial (final) state. The functions and /2(w) are fixed in terms of the 



dark-matter electric charge and magnetic moment, respectively, as w ^ [U, |19|, ioj with- 
out further assumption — it does not even matter if the dark-matter particle is composite. 
The amplitude A4r{k,p — k,p) is implicitly a 2 x 2 matrix in the photon polarization, and 
its diagonal matrix elements describe dispersion in propagation and attenuation [l^. The 
f2{(^) term describes changes in polarization with propagation, so that we need not consider 
it further. Under analyticity and unitarity, expanding fi{u) for u ^ Uth yields a series in 
positive powers of u"^ for which the coefficient of every term of 0(c<;^) and higher is positive 
definite [ll|, Il2| . Thus a term in n{uj) which is linear in co, discussed as a signature of Lorentz 
violation [9|, l21j, does not appear if < coth and the medium is unpolarized. We parametrize 
the forward Compton amplitude as Air = ^j=o^2j^'^'' , where Aq = —2e^e^ [ij, [l7(] and 
the dark-matter millicharge is ee. The terms in OiuJ^) and higher are associated with the 
polarizabilities of the dark- matter candidate. 

Dispersive effects in light propagation are controlled by the group velocity fg, 
light emitted from a source a distance / away has an arrival time of tiiS) = }ly. 
distant sources we must also take the cosmological expansion into account |22|], so that as 
we look back to a light source at redshift z, we note that the dark-matter density accrues a 
scale factor of (1 + 2;)^, whereas the photon energy is blue shifted by a factor of 1 + z relative 
to its present-day value uq Thus the light arrival time ^(cjo, z) is 
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with the Hubble rate H{z') = Hq^J (1 + z'^VLm + ^a- We employ the cosmological param- 
eters determined through the combined analysis of WMAP five-year data in the ACDM 
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model with distance measurements from Type la supernovae (SN) and with baryon acous- 
tic oscillation information from the distribution of galaxies jsj. Thus the Hubble constant 
today is Hq = 70.5 ± 1.3 km s~^Mpc~^, whereas the fraction of the energy density in matter 
relative to the critical density today is Qm = 0.274 ± 0.015 and the corresponding fraction 
of the energy density in the cosmological constant A is = 0.726 ± 0.015 [3|. 

We find that the time delay is characterized by powers of u"^ and unknown coefficients A2j. 
Different strategies must be employed to determine them. The A2 term incurs no frequency- 
dependent shift in the speed of light, so that to infer its presence one needs a distance measure 
independent of z, much as in the manner one infers a nonzero cosmological constant from 
Type la supernovae data. Interestingly, as > it has the same phenomenological effect as 
a nonzero cosmological constant; the longer arrival time leads to an inferred larger distance 
scale. Cosmologically, though, its effect is very different as it scales with the dark-matter 
density; it acts as grey dust. The remaining terms can be constrained by comparing arrival 
times for differing observed Uq. 

The determination of Aq and A4 require the analysis of the GRB light curves at extremely 
low and high energies, respectively, and probe disjoint dark-matter models. As a first ap- 
plication of our method, we use radio afterglow data to determine Aq and thus to yield a 
direct limit on the electric-charge to mass of dark matter. This quantity gives insight into 
the mechanism of dark matter stability. If dark matter possesses an internal symmetry, e.g., 
it cannot decay to lighter particles and conserve its hidden charge. Such dynamics can also 



conspire to give dark matter a slight electric charge |23H25| . which, no matter how small, 
reveals the existence of its hidden interactions and the reason for its stability. 

To determine Aq we consider GRBs with known redshift in which a radio afterglow is 
also detected. We collect the data and describe the criteria used in its selection in the 



supplementary material [26|. The time lag between the initial detection of the GRB at 
some energy and the detection of the radio afterglow is r = t{uj^Q'^,z) — t{uQ^^^,z). If we 
first observe the GRB at keV energies and compare with the observed arrival time in radio 
frequencies, then the terms in positive powers of Uq, as well as the term in 1/(cjq'^'^)^, 
are negligible; we let In order to assess reliable limits on Aq we must separate 

propagation effects from intrinsic source effects. Statistically, we expect time delays intrinsic 
to the source to be independent of z, and the time delay from propagation to depend on z 
and w in a definite way Such notions have been previously employed in searches for Lorentz 
invariance violation j21|. We separate propagation and emission effects, respectively, via 

= Ao^ + 6iil + z)u) (2) 

1 + z 

where K(z) = (1 + zY^ dz (1 + z')H{z')~^ depends on the cosmological past through the 
Hubble rate H{z) and 5((1 + z)!/) allows for a frequency-dependent time lag for emission 
from the GRB in the GRB rest frame. The frequency u = u/27c, and Ao contains the 
millicharge-to-mass ratio e/M, i.e., Att'^Aq = —Aopo/AM'^ = 2Txae'^ po/ M"^ with po — 1-19 x 
10~^ GeV/cm^ |3i] and a the fine-structure constant. To provide a context, we first consider 
the value of |e|/M which would result were we to attribute the time lag associated with the 
radio afterglow of one GRB to a propagation effect. Choosing the GRB with the largest 
value of K{z)/i''^, we have a time lag of 2.700 ± 0.006 day associated with GRB 980703A at 
z = 0.967 ± 0.001 measured at a frequency of 1.43 GHz. With Eq. ([2]), setting 6 = 0, and 
noting that K{z)/p^ = 1170 ± lOMpc GHz'^ if the errors in its inputs are uncorrelated, 
the measured time lag fixes l^l/M ^ 9 x 10^^ eV^^. Since there are no known examples of 
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FIG. 1: The time lag r determined from the observation of a GRB and its radio afterglow plotted 
as a function of K{z)/u'^ with u = lo /27r, employing the data reported in the supplementary 
material 
GHz (I 
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The points correspond to frequency windows of 4.0 — 12 GHz (t, green), 12 — 30 
maroon), and 30 — 75 GHz blue) in the GRB rest frame. Points with (1 + z)iy < 4.0 
GHz do not appear within the chosen frame of the figure. The fit of Eq. ([2]) to the data with 
(1 + z)!^ > 4.0 GHz with a scale factor in the uncertainty in t/(1 + z) of 450, to compensate for 
environmental effects in the vicinity of the emission from the GRB, yields Aq = O.OOlO it 0.0019 
day Ghz2 Mpc'^ and 6 = 0.65 ± 0.10 day with x^/ndf =1.13. Thus Aq < 0.005 day Ghz^ Mpc^^ 
at 95% CL to yield |e|/M < 1 x 10~^ eV~^ at 95% CL. The statistical scale factors are not shown 
explicitly. For clarity of presentation we display time lags in the GRB rest frame of less than 10 
days only. 



a radio afterglow preceding a GRB, this single time lag in itself represents a conservative 
limit. Turning to our data sample of 53 GRBs, we plot the measured time lag versus 
K{z)/v'^ in Fig. [Hand make a least-squares fit of Eq. (j2]) to determine and + z)v). 
We require Aq > Q as demanded by our model. Fitting to the points with frequencies of 
4.0 - 75 GHz in the GRB rest frame, we determine \e\/M < 1 x 10"^ eV"^ at 95% CL, 
which is comparable to our limit derived from a single observation of GRB 980703A. The 
dependence of our fit results on the selected frequency window, as well as the stability of our 
fits to the significance of the radio afterglow observation, to evolution effects in z , and to the 
more poorly determined red shifts and radio afterglows is discussed in the supplementary 
information (26| . 

We have found a direct observational limit on the dark-matter electric-charge-to-mass 
ratio. Our study probes for a charge imbalance averaged over cosmological distance scales, 
without regard to its sign, at distance scales shorter than the wavelengths of the radio 
observations in our data set. Our bound rules out the possibility of charged "Q-balls" |25l. l27 



of less than 100 keV in mass as dark-matter candidates. Our limit holds regardless of the 
manner in which the dark matter is produced, though we can compare it to limits arising 
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from the nonobservation of the effects of milhcharged particle production. For example, 
for M ~ 0.05 eV, they are crudely comparable to the strongest bound from laboratory 
experiments for |e| < 3 — 4 X 10 ^ for M < 0.05 eV |29|]. In comparison the model- 

independent bound arising from induced distortions in the cosmic microwave background 
(CMB) radiation is \e\ < 10^^ for M < 0.1 eV, though model-dependent constraints reach 



e\ < 10 ^ for M < 2 X 10 ^eV 30|. Cosmological limits also arise from observations of 



the Sunyaev-Zeldovich effect for which, e.g., \e\ < 3 x 10~^ for M ~ 10~^eV [31]. For 
these light masses our limit is stronger, which shows that millicharged particles of such mass 
and charge are not the primary constituents of dark matter. Limits also arise from stellar 
evolution and big-bang nucleosynthesis constraints, for which the strongest is |£:| < 2 x 10~^^ 
for M < 5keV [28|, as well as from the manner in which numerical simulations of galactic 



structure confront observations 
Fig. [2J Indirect limits can be evaded: 



We offer a visual summary of this discussion in 
for example, in some models, the dynamics which 

other 



gives rise to millicharged matter are not operative at stellar temperatures [29|, |35 
models evade the galactic structure constraints (25I . 27 1. 

Our limit also significantly restricts the phase space of models with hyp erweak gauge inter- 
actions and millicharged particles which can arise in string theory scenarios jsi, Sgj as viable 
dark-matter candidates. We estimate that our limit can be improved considerably before the 
dispersive effects from ordinary charged matter become appreciable The largest such 
contribution to Aq should come from free electrons. We estimate the cosmological free elec- 
tron energy density pe to be no larger than = (Me/Mp)pcr^6 ~ 0.130 eV/cm^ [3j, where ilb 
is the fraction of the energy density in baryons with respect to the critical density today and 
Me and Mp are the electron and proton mass, respectively. Replacing po with pe and e/M 
with l/Mg in Aq we find that our limit would have to improve by 0{2 x 10^'^) before the con- 
tribution from free electrons could be apparent. We set our limit of |e|/M < 1 x 10 ~^ eV~^ 
from existing radio observations at no less than 4 GHz in the GRB rest frame, so that our 
limit is certainly operative if ci;th/2vr > 4 GHz, or crudely, if M > 8 x 10~^ eV. Studies of 
the polarizability in QED for which Uth = 0, also reveal the analytic structure in uj 
we have assumed for the forward Compton amplitude. Thus we believe our limit to be of 
broader validity, so that the lower limit on the mass can be less than 8 x 10~^ eV, though 
it is model dependent and set by the Lee- Weinberg constraint 
mass of ~ 6 X 10~^ eV is determined in the axion model of Ref. 
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much as the minimum 
. Forward scattering is 

coherent irrespective of whether the photon wavelength is large compared to the interpar- 
ticle spacing, so that we expect our results to persist in the dilute particle limit as well, as 
supported by laboratory studies 13]. One further comment: at a frequency of 4 GHz our 
limit implies that we probe the average net charge of dark matter, with no constraint on 
its sign, at length scales of no longer than 8 cm. Our limits can be significantly bettered 
through GRB radio afterglow studies at longer wavelengths. 

We thank Keith Olive for an inspiring question and Scott Dodelson, Renee Fatemi, Wolf- 
gang Korsch, and Tom Troland for helpful comments. SG would also like to thank Stan 
Brodsky for imparting an appreciation of the low-energy theorems in Compton scattering 
and the Institute for Nuclear Theory and the Center for Particle Astrophysics and Theo- 
retical Physics at Fermilab for gracious hospitality. This work is supported, in part, by the 
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FIG. 2: Comparison of our direct limit on the absolute electric charge of dark matter, \e\, in units 
of e with candidate mass M with limits from other sources stemming from millicharged particle 
production. Our limit derives from observations of GRBs and their radio afterglows and is marked 
"ORB" (solid, black), so that the region above that line is excluded. The condition u < ujth, as 
discussed in text, sets the lower endpoint of the solid line. Since Eq. ([2|) persists in QED, for 
which ujiii = 0, we expect our limit to persist for lighter, cold dark matter as well, as indicated 
by the dashed line. The strongest laboratory limits, which are for fermions, are marked "Pol" 
(solid, blue) 29|, and the strongest limits from induced distortions in the CMB, which are also for 
fermions, are marked "CMB" (long-dashed, maroon) — the upper curve is the model-independent 
limit, whereas the lower curve is the model-dependent (md) limit [30|. Constraints on \£\ emerge 
from limits on novel energy-loss mechanisms in stars and supernovae; such limits also fail to act 
if |ej is too large. The limit from plasmon decay in red-giants is marked by "RG" (short-dashed, 
red) [2S], the same limit in white dwarfs is marked by "WD" (dot-dot-dashed, indigo) [2^, and 
the limit from SN 1987A is marked by "SN" (dot-dashed, orange) |2a]. The RG limit acts if 



^1 ^ 10 ^ 28|- We have also reported the limit from big-bang nucleosynthesis, marked by "BBN" 



(dotted, green), from Ref. [28|] as well 



U.S. Department of Energy under contract DE-FG02-96ER40989. 
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Supplementary Information 



I. GAMMA-RAY BURST DATA 

From the publicly available data, we consider the gamma-ray bursts (GRB) with known 
redshifts and detected radio afterglows through March, 2009. Tables 1-4 collect the GRBs 
with these properties which satisfy certain criteria. We demand that the energy of the initially 
detected gamma-ray burst in its rest frame be compatible with the energy range of the 
Fermi Gamma-Burst Monitor. This excludes the X-ray flashes GRB 080109A and GRB 
020903 A. To determine whether an observation of the radio flux is indeed a detected radio 
afterglow, we demand that the radio flux detection be within the error box of the location 
of the observed GRB, that it not be a site of previously observed radio activity, and that 



it be significantly non-zero. We use the criterion employed by Ref. [55|] and thus require 
that the radio observation be nonzero by three standard deviations or more to be termed a 
detection. Such considerations exclude GRB 030277A, GRB 980425A, and GRB 011130A. 
Consequently, we find 53 GRBs to consider, for which we report all detected radio frequencies 
of 75 GHz or less. We report the time of the initial detection of the GRB and the time of the 
detection of a particular observation frequency v in the associated radio afterglow (RA). We 
calculate the time lag between the detection of a particular radio frequency and the initial 
detection of the GRB from these quantities. 

Remarks concerning the collected data are in order. A redshift marked with an asterisk * 
means that no certainty has been reported in the work which determines it. In those cases, 
we assume the uncertainty in the redshift to be plus or minus one unit in the last significant 
figure reported . The redshift we report for GRB 021004A is the average of those reported in 
Refs. 



145 



1461 ]. For the redshifts marked with a circle °, the associated publication states 
that the reported redshift is either a lower limit or the "likely" value. In the particular case 
of GRB 970 508A, work in refinement of the upper bound on the red shift has been reported 
in Ref. 204]. For the frequency marked with a double dagger -f, observations were made in 
the frequency range from 14.5 GHz to 17.5 GHz. 

We also note which radio detections have lesser significance, so that detections with a signif- 
icance between 3(j and 4(T are marked with a sharp ^, whereas detections with a significance 
between 4f7 and 5(T are marked with a flat ^ For time lags marked with a dagger we assume 
the uncertainty in the time lag to be plus or minus one unit in the last significant figure of 
the reported detection time of a particular radio frequency. If a radio frequency detection 
reports an observation time interval, the time lag is calculated from the midpoint of the time 
interval, and we take one half of the observation time as the uncertainty in the time lag. The 
RA time reported in these cases is that of the midpoint of the observation time. 
To realize the limits reported in the main body of our paper, we fit to Eq. (2) of that source. 
We find that large statistical scale factors must be employed to realize good fits to the 
data. This may stem, in part, from the circumburst environment, and, more generally, from 
time delay effects which arise from neither source effects nor propagation effects across the 
expanse of s pace. The nature of the circumburst environment of GRB 050904A is discussed 



in Ref. [111]. 



To study how our limits rely on the details of our data set, we compare fits with frequency in 
the source rest frame. We also study how the fits change with the significance of the included 
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radio observations, with omitting radio observations with reported integration times, and 
with omitting more poorly determined red shifts. We also try to study source evolution 
effects with z. As to the frequency dependence, differences do emerge if we include points at 
the observed frequency of 1.43 GHz. This portion of the data set admits time delays in the 
GRB rest frame of a month and more, and larger scale factors are required to yield fits of 
comparable quality. Fitting data satisfying (1 + z)v > 4.0 GHz, for which with a scale factor 
of 450 applied to the error in r/(l + z) we find xV^df = 1.13 with Aq = 0.0010 ± 0.002 
day Ghz2 Mpc"^ and 5{{1 + z)u) = 0.65 ± 0.10 day, or Aq < 0.005 day Ghz^ Mpc'^ at 95% 
CL. Were we to fit all the data points we would require a scale factor of 685 to yield a fit of 
X^/ndf = 1.13 from which we would find Aq < 0.007 day Ghz^ Mpc"^ at 95% CL. We thus 
restrict our discussion henceforth to the portion of the data set which satisfies (1 + z)z^ > 4.0 
GHz. If we now restrict our fit to points in which the significance of the radio observation 
is 5(7 or more, then our limit on the slope is still Aq < 0.005 day Ghz^ Mpc^^ at 95% CL, 
though we must increase the scale factor to 512 to recover a fit of comparable quality. Some of 
the RA measurements have reported integration times and hence determinable uncertainties, 
whereas others do not. Repeating the fits without the points with reported integration times 
yields no significant difference in the fit, though we must increase the scale factor to 467 to 
yield a fit of comparable quality. This may well follow as the dropped points do have larger 
errors and hence play a lesser role in the original fit. If we now return to our original data set 
and omit those points for which the determined redshift is a lower limit or merely "likely," 
the limits do weaken. For example, with a scale factor of 419, we find x^/ndf = 1.13 with 
Ao = 0.0053 ± 0.0048 day Ghz^ Mpc'^ and 6{{l + z)u) = 0.41 ± 0.21, so that the limit 
on the slope at 95% CL is roughly a factor of 3 worse. This can be rationalized as these 
poorly known redshifts possess large z. Finally, to study the effects of cosmological evolution 
on the GRBs, we repeat our fit with z > 1 and find similar results. That is, for a scale 
factor of 539 we determine x^/ndf = 1.14 with Aq = 0.0017 ± 0.0023 day Ghz^ Mpc'^ and 
5{{1 + z)u) = 0.58 ± 0.12, to yield Aq < 0.006 day Ghz^ Mpc'^ at 95% CL. Thus the effects 
of cosmological evolution appear to be small. 
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19:55:54.480 [107] 
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TABLE I: Gamma ray burst and radio afterglow data from mid-2005 through March, 2009. 
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TABLE II: Gamma ray burst and radio afterglow data from mid-2002 through mid-2005. 
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GRB 


Redshift 


GRB Time 


RA Time 


Frequency 


Time Lag 






(UT) 


(UT) 


(GHz) 


(d) 


020405A 


0.695 ± 0.005 [155] 


00:41:26 [156] 


Apr 06.22 [157] 
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1.19 ± O.Olt 


020405A 






Apr 08.38 \U] 


4.86" 
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2.142 ± 0.001* [158] 


19:09:21 [159] 
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22.5" 


7.76 ± O.Olt 
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0.362 ± 0.001 [161] 
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Nov 22.83 [163] 


8.70 


1.05 ± O.Olt 
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Nov 25.20 [163] 
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3.42 ± O.Olt 


010921A 


0.450 ± 0.005 [164] 
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Oct 17.15 [166] 
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25.93 ± O.Olt 
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25.93 ± O.Olt 


010222A 


1.477 ± 0.001* [167] 


07:23:30 [168] 


Feb 22.62 [169] 


r 

22'' 


0.31 ± O.Olt 


010222A 






Feb 23.66 [170] 


8.46 


1.35 ± O.Olt 


010222A 






Feb 24.56 [170] 


4.861* 


2.25 ± O.Olt 


000926A 


2.066 ± 0.001*° [171] 


23:49:33 [172] 


Sep 28.17 [173] 


8.46 


1.18 ± O.Olt 


000926A 






Sep 29.726 [173] 


4.86 


2.733 ± O.OOlt 


000926A 






Oct 04.186 [173] 


22.5 


7.193 ± O.OOlt 


000911 A 


1.0585 ± 0.0001 [174] 


07:15:25 [175] 


Sep 14.36 [174] 


8.46 


3.06 ± O.Olt 


000911 A 






Sep 22.35 [174] 


4.8611 


11.05 ± O.Olt 


000418A 


1.1181 ± 0.0001 [176] 


09:53:10 [177] 


Apr 29.07 [178] 


8.46 


10.66 ± O.Olt 


000418A 






May 03.04 [178] 


4.86 


14.63 ± O.Olt 


000418A 






May 03.04 [178] 


22.46 


14.63 ± O.Olt 


000418A 






May 07.18 [178] 


8.35 


18.77 ± O.Olt 


000301C 


2.0335 ± 0.0003 [179] 


09:51:37 [180] 


Mar 04.98 [181] 


15.0^ 


3.57 ± O.Olt 


000301C 






Mar 05.67 [181] 


r 

4.86'' 


4.26 ± O.Olt 


000301C 






Mar 05.67 [181] 


8.46 


4.26 ± O.Olt 


000301C 






Mar 05.67 [181] 


22.5'' 


4.26 ± O.Olt 


991216A 


1.00 ± 0.02 [182] 


16:07:01 [183] 


Dec 17.783 [184J 


4.8 


1.11 ± 0.13 


991216A 






Dec 18.00 [185] 


15^ 


1.33 ± O.Olt 


991216A 






Dec 18.16 [185] 


8.46 


1.49 ± O.Olt 


991216A 






Dec 18.32 [185] 


8.42 


1.65 ± O.Olt 


991216A 






Jan 03.11 [1851 


4.86'' 


17.44 ± O.Olt 



TABLE III: Gamma ray burst and radio afterglow data from mid-1999 through mid-2002. 
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TABLE IV: Gamma ray burst and radio afterglow data from 1997 through mid-1999. 
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